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Abstract. In this note, we compute the Hadamard coefficients of (algebraically) 
integrable Schrodinger operators in two dimensions. These operators first ap- 
peared in |BL| and [B] in connection with Huygens' principle, and our result 
completes, in a sense, the investigation initiated in those papers. 



Let C, = — A„ + be a Schrodinger operator on R" (n > 1) with C°°-sniooth 
potential V — V{x) defined in some open domain C M". Recall that the heat 
kernel of C is the solution <i>_|. {x, ^, t) e C°° (flx^lx R^) of the initial value problem 

(1) (^^^+£^^+{x,^,t)^0 , \un^^+{x,^,t)=S{x~0 , 

where 6{x — is the Dirac delta-function on Q with support at ^. It is well known 
(see, e.g., [R], Sect. 3.2.1) that <i>+(a;,^,t) has an asymptotic expansion of the form 

(2) <i>+ix,tt) - [l + J2uAx^0t''j ast->0+, 

with coefficients U^{x, ^) e C°°{Q x f2) determined by the following transport equa- 
tions 

(3) ix-C,d,)U,ix,0+iyU,ix,0^-^U,-ii-:mx) , iy = l,2,... 

The system ^ has a unique solution {[/^/(a;, ^)}j^Q if one sets Uo{x,£^) = 1 and 
requires each Ui^(x, £,) to be bounded in a neighborhood of the diagonal x = ^. Fol- 
lowing [G], we will refer to {U^{x, ^)} as the Hadamard coefficients of the operator 
C. 

In general, calculating the Hadamard coefficients for a given potential V is 
a difficult problem. Of special interest are potentials, for which the heat kernel 
expansion ^ is finite, i. c. the sum in the right-hand side of ^ has only finitely 
many nonzero terms. In this case, formula ^ yields not only a "short-time" 
asymptotics, but an exact analytic representation for ^+{x, t) valid for all t G M.^. 
Such potentials are usually called Huygens ' potentials in view of an important role 
they play in the theory of Huygens' principle (see, e. g., [B Vj . Ch. I). 

The problem of describing all Huygens potentials goes back to Hadamard's clas- 
sical treatise jHl and still remains open in all dimensions, except for n — 1 (see 
[L]) and n = 2 (see [B ). In dimension one, these potentials coincide with the 
well-known Adler-Moser potentials |AMj . which are the rational solutions of the 
Korteweg-de Vries hierarchy of nonlinear integrable PDE's (see [S] or [By|, Ch. 3, 
Sect. 3.2); apart from the original works [LSj and [Lj, their Hadamard coefficients 
have been studied recently in [Grj , [11] , [E] and [Ha| . 



Berest's work partially supported by NSF grant DMS 04-07502. 

1 



2 



YURI BEREST, TIM CRAMER, AND FARKHOD ESHMATOV 



In the present paper, we will deal with the two-dimensional Huygens potentials, 
which have received so far much less attention (see, however, [CFV| ). Our main 
result (Theorem [5]) provides simple and explicit formulas for the Hadamard co- 
efficients of the corresponding Schrodinger operators. It is surprising that these 
formulas do not seem to have analogues in dimension one: in a sense, the two- 
dimensional case is simpler than the one-dimensional one! 

We begin by recalling the main result of BL^ and [B]. Let k = {ko, ki, ... ,km) 
be a finite, strictly increasing sequence of integers, with /sg = , and let ipi be real 
numbers given one for each integer ki , with tpo = . Passing to the polar coordi- 
nates {xi, X2) = (rcosi/J, rsiuip) in K^, we associate to these data the following 
potential 



(4) Vk{xi,X2) = — - — logWr[xo, xi, ■ ■ ■ , Xr. 



2 

where Xi ■= cos{kiip + (pi) , i = 0, 1, ...,m, and Wr[xo, Xi: •■■ . Xm] is the 
Wronskian of the set {xi} taken with respect to the variable (p. As all fc^'s are 
integers, 14 is a single- valued rational function on R'^, homogeneous of degree 
—2, whose analytic continuation to has singularities along certain lines passing 
through the origin. For example, for k = (0, 1, 3, 4) , with all (/s^'s being 0, we have 

UiA9xt + 28xlxl-x*) 

Vk[Xi,X2) - 



4 {^A + A? 

In general, ([4]) depends on both the choice of kis and the choice of (^^'s, though 
we suppressed the latter from our notation. 

Theorem 1 ( [BU lB ] ). A (locally) smooth function V onM.'^, which is homogeneous 
of degree —2, is a Huygens potential if and only if V = Vk{xi, X2) for some integer 
sequence k = (ki) and real numbers {(fi). 

Remark. The "if" part of Theorem [T] was first proven in BL and then reproven 
by a different method, together with the "only if" part, in [B] (see loc. cit., The- 
orem 1.1). The assumption that V is homogeneous can be relaxed and replaced 
by a weaker condition that V and all the Hadamard coefficients of V are algebraic 
functions (see [CFV| ). 

The main result of this paper can now be encapsulated in 

Theorem 2. The Hadamard coefficients of C = — A2 -f Vk with potential ^ are 
given, in terms of polar coordinates x = {rcosip, rsiiup) and ^ = {g cos cf), gsincf)) , 
by 



(5) c/,(x,0 = ^^E^^*''(^)*''('^)^ir^(cos(^-'^)) ' Vz.>0, 
where 



{rgY ^ 



(6) Wr[x.,x.,.-..X.-uX.^u-.-,Xra] ,^ , _ 2) 

Wr[xo, Xi, ■■■ fj-g 

andTniz) := cos(iV arccos z) is theN-th Chebyshev polynomial, with t\^\z) being 
its derivative of order v with respect to z . 
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Remark. It follows immediately from ^ that U,y{x,£,) = for i/ > km , implying 
that Vfc is a Huygens potential. For k of length one, i. e. k — (0, N) , formula 
(H]) yields the Calogero-Moser potential of dihedral type l2{N), and in this special 
case the coefficients ([5]) have already been found in |BL] (see he. cit., Sect. IV, 
Example 1). 

There are several ways to prove Theorem [2] Perhaps, the most straightforward 
one is to use a differential recurrence relation between the Hadamard coefficients 
of operators Ck and C-^ , where k — [k, fcm+i) is obtained by adding one integer 
on top of k (see [BL' , (82)). This method requires double induction (in m and v) 
and leads to rather unwieldy calculations. 

Here, we will offer a more illuminating argument based on the remarkable fact 
that the coefficients Ui,{x,S,) appear not only in fundamental solutions of the heat 
equation but also in its elliptic and hyperbolic counterpart^^- Instead of the Cauchy 
problem ([T|), we will consider 

(7) /:[G{-,m^)^s{x-o , Gi-^o^-D'm, 

where I>'(f7) denotes the space of distributions on C°°-functions with compact 
support in C R" . Of course, unlike the heat kernel, G{x, ^) is not uniquely 
determined by ([7]), but only up to adding smooth functions from Ker(jC) . The 
problem is now to describe the singularities of G{x,^) . In modern language, the 
solution to this classical problem is given in terms of Riesz distributions, and as 
in the hyperbolic case (see, e.g., [D]), it depends on whether n is even or odd. 
Specifically, if 71 > 2 is even, G{x, £,) has the following asymptotics "in smoothness" 
(cf. [Hg, Sect. 3, (1.4)): 

00 

(8) G(x,0-^E^-(^'^)^^-^(t) ' 

where 7 ~ \x — ^\'^ is the square of the Euclidean distance between x and ^, 
S\{t) := t^/r{X + 1) is the family of Ricmann-Liouville distributions on R^, de- 
pending analytically on the parameter A G C, and 5*^^ := dS\/dX are the adjoint 
distributions of S\ with respect to A (see |GSj . Ch. VI, Sect. 2). The coefficients 
Ui,{x,£^) in ([H are smooth functions, which satisfy, up to rescaling factor —1/4, 
the same transport equations ([3]) as the heat kernel coefficients U^{x,S^) . Like in 
the heat kernel case, this can be verified by direct calculation, substituting ([5]) into 
([7]). By uniqueness of the regular solution of ([3]), wc thus have 

(9) U,ix,0^ (-7] U^^^O foralli.>0. 



Now, let n = 2 and assume that V — V{xi, X2) is (locally) analytic, as are 
our potentials (|4]). In this case, the distributions S'^ appear in ([8]) only with non- 
negative integer A's, and for such A's they can be easily calculated 

(10) Slit) : 



dX 



= ^f1\ogt + C^f1 , t^^O, 1, 2, 



where Ci, G R are some constants, with Co = 0. Substituting ^ and ITUl) in 
we get 

G{x,0-W{x,0 log7+ ••■ , 



^For an excellent survey on this classical subject we refer the reader to |Ba| . 
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where " . . . " denote smooth functions, which do not contribute to singularities of 

G{x,S,), and 

oc 

(11) W^(^-0:-E7ZivV^-(^'^)^''- 

Since V is analytic, the series uniformly converges in a neighborhood of the 
diagonal x = ^. This was already observed by Hadamard, who called W{x,^) the 
logarithmic term of the "elementary solution" G{x, ^). He also showed that W{x,^) 
is the (unique) analytic solution of the Goursat problem 

(12) C[W{-,i)]{x)=0 , W\^^o^l, 

with the boundary condition given on the (complex) characteristics of £ (see [H] , 
Ch. Ill, Sect. 46). 

Now, to prove Theorem [2] we will simply solve (fT2l) for the operator C with 
potential ([4]), using the method of separation of variables, and then recover the 
coefficients Ui,{x,i^) by expanding the solution in the vicinity of 7 = as in (|lip . 

First, we introduce a one-dimensional Schrodinger operator L by writing C in 
terms of polar coordinates (cf. [B], Sect. 3): 

Explicitly, L = —d^jd^P' + 14 (cos 1^9, sin(^) , where 14 is given by (jj]). Next, we 
prove 

Lemma 1. 7/5'== ^'('p) satisfies L['i>] = fc^^ with some integer k > 0, then 
'i'{(p)Tk[{r/g + g/r)/2] e KeT{C) for any ^ . 

Proof. Changing the variables (r, ip) (z, tp) in (|13p . with z := (r j g -\- Q/r)/2 , 
yields 

If F)"^ rl ' 

^=4 (i--')#2--l^+^ 

Since Tk{z) is an eigenfunction of [1 — z'^)d?/dz'^ — zd/dz with eigenvalue —fc^, the 
claim is obvious. □ 

Lemma 2. The functions '^i defined by ([6]) satisfy 

(14) L[^>,] - fc2 , I = 0, 1, . . . , m , 

(15) Ec,;^',((p) *,((/)) COs{h{ip - (f))) = 1 . 



Proof. The equations (fTi|) simply say that 5'i's are eigenfunctions of L with 
eigenvalues This follows directly from Crum's classical theorem (see [C], p. 124). 

The identity (jlSp seems more interesting: we could not find it in the literature, 
so we will prove it in detail by induction on m. 

The case m = 1 is straightforward. Writing S'm(iy9, </>) for the left-hand side of 
(lisp , we now fix m > 1 and assume that Sm{'^, 0) = 1 for all sequences of integers 
(fco = 0, fci, . . . , fc„i) and reals ((/Jq = 0, ipi, . . . , Lp^) of length m . To make the 
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induction step we add km+i G Z , km+i > km , and Lpm+i £ M to these sequences 
and consider 



TJJ.+ l 



i=0 

where and Ci are defined by formulas ^ with m replaced by rn + 1 . By Crum's 
Theorem, ^I^i's are eigenfunctions of an operator L obtained from L by applying 
the Darboux transformation (cf. "B], (3.29)): 

(17) L ~ A*^ O Am + km+l ^ L = Am O A*„ + fc^_|_i , 



where 



(18) Am := *,„Vi o i^—j o ^m+i , v4:„ := o i^ — j o ^-^^i ■ 

It is immediate from (fT8|) that . On the other hand, we have 

(19) = for aU i = 0, . . . , m . 
In fact, a trivial calculation shows that (fT9|) is equivalent to 

d f W, \ WW,,m+l 



(20) 



\ Wm+lJ ~ 



m+1 



where W denotes the Wronskian of the set {xo, Xi, ■ • ■ , Xm+i}, and Wi , Wi^m+i are 
the Wronskians of this set with functions Xi and {xi, Xm+i} being omitted. Now, 
to prove ((20|) consider the system of linear equations 



(21) E^»^'V» = xL'li , ^ = 0,1, 

i=0 



By Cramer's Rule, the solution to this system is given by j/i = (— 1)™+» Wi/Wm+i ■ 
On the other hand, differentiating both sides of l|2T|) with respect to (f yields 

E"oxf^2/- = for fc = 0, 1, m - 1 and E" o X^^ J/- = W/Wm+i for 
k = m. Solving these equations for ?/■ , we get y[ — {-'l)™^'^WWi^m+i/Wm+i , 
which is equivalent to ([20|) . 

It follows from that ^^['1'^] = -~AmA*^[^i] = ~{L - fc^^i)*^ . Hence 

(22) A^J^-,] = (fc^+i - fcf)^', for alU = 0, . . . , m . 

Now, applying Am and yl*^ to p6|) . we can formally relate 5to+i(93, 0) to Smi'P, 0)- 
In fact, a straightforward calculation using and (|^^ shows 

(23) (•,(/.)] + [5,„+i iip,-)]^Am[Smi-, <(>)]+ Am [Sm (V,-)], 

where omitted are the variables on which the differential operators act. Since 
Smi'P, 0) = 1 by induction assumption, we can regard as a first order PDE for 
the function 5m+i(</5, 0)- Specifically, substituting Sm — 1 into (j23p yields 

f d d \ ^ fdB dB\ „ dB dB 

where B := log m+iif)'^ m+i{4')] ■ This PDE can be easily integrated: changing 
the variables x := {ip + (f)) / 2 , t [tp — (j))/2 and F := {Sm+i — 1) transforms 
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([24|) to the equation dF/dx = , which means that F depends only on t. It follows 
that 

(25) S^+i {lp,(I))^1 + F{lp- 4>) ^,n+i{v) ^,n+l{4>) , 

where is a differentiablc function of one variable defined on (0, 27r). 

Now, if we replace km with km+i in the sum Sm and repeat the above argument, 
adding km to the partition (ko = 0, ki, . . . , km-i, fcm+i) , then, instead of we 
get 

(26) Sm+ii^,cb) - l + G((^-(/.)^„((/j)^„((/.) , 

where G is another differentiablc function on (0, 27r). Comparing ((25)l and ((26|) 
shows that F — G = 0. Indeed, if one of these functions {F say) is nonzero, then 
by continuity F{ip — (f)) 7^ in an open subset of (0, 27r) x (0, 2it) , and in that 
subset we have 

(27) G((^-0)/F(^-0) =/i(^)/i(0) , 

where h :— ^m+i/'^m ■ Differentiating both sides of (P7|) with respect to ip and (f) 
and adding the results yields h'{ip) h{(f)) + h{Lp) h'{(j)) = . Whence, letting = </> , 
we see that h must be constant. This means that V&m+i is a multiple of V&m, which 
is impossible, since and are nonzero eigenf unctions of L corresponding 

to different eigenvalues (fc^ and k^+i respectively). Thus F = 0, and therefore 
Sm+ii'^T (t>) = 1, finishing the induction. □ 

Now, combining the results of Lemmas [T] and [21 we see at once that 

(28) ly c,*,(^)vE',(0)rfeJ(r/f?+g/r)/2] 

i=0 

is a solution to ([72]) . Indeed, by ([14]) and Lemma [1] each summand of ([28)1 lies in 
the kernel of C, and hence C[W] = by linearity of £. On the other hand, in polar 
coordinates 7 = + — 2 rg cos{ip — (j)) , so 

(29) U- + ^)=^+co8{^-cj)), 
2 \g r J 2rg 

and therefore 

TkA{r/g + g/r)/2] ^ cos{h{^~(j))) on 7 

for all i — 0, 1, ... , m . The boundary condition for W follows then from ([TS]) . 

Now, by [B , Lemma 3.1, U^{x,S,) are homogeneous functions of x and ^, which 
can be written in terms of the polar coordinates as 

(30) U,ix,0 = -^MV,^) , i'>0. 

Clearly, there is at most one expansion of W{x,^) of the form ([TT|) with coefficients 
([30|l . To find this expansion, we take the obvious Taylor formulas (see ([29| ) 

(31) T,A{r/g+g/r)/2] = J2- (^) Tif)(cos(^ - 0)) , 

substitute (|3ip into ([^5)1 and reorder summations. Collecting coefficients under the 
different powers of 7 gives then the desired formulas ([5]). 
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In the end, we note that our operators C are examples of algebraicahy integrable 
Schrodinger operators (see [CV| ). Such Schrodinger operators possess special eigen- 
functions called the Baker- Akhiezer functions. Like fundamental solutions above, 
the Baker- Akhiezer functions have asymptotic expansions, with coefficients sat- 
isfying (up to rescaling factor 1/2) the same transport equations ([3]). This was 
originally discovered in the special case of Calogero-Moser potentials in [BVlj , but 
the argument of |BVlj applies to any homogeneous operator (see, e.g., [CFV] ) . In 
combination with Theorem [51 this yields 

Theorem 3. The Baker- Akhiezer function of C = — A2 + Vk{xi^X2) is given by 

where Ui,{x, ^) are the same coefficients as in ([5]). 
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